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摘 要 :高 光谱 图 像 (Hyper Spectral Image，HSI) 分 类 是 HSI 人 处理 中 的 重要 预 处 理 手段 ,其 目标 是 对 HSI 数 据 中 每 个 像素 
点 进行 类 别 标 记 , 标 记 结 果 常 用 于 识别 .勘探 等 应 用 。 针 对 HSI 分 类 任务 中 存在 的 数据 量 大 .数据 维度 高 .已 知 样本 量 少 等 
难点 ,提出 一 种 基于 图 模型 的 半 监 督 分 类 算法 。 该 算法 将 HSI 数 据 建 立 为 图 以 实现 降 维 ,而 后 将 分 类 问题 归结 为 一 个 无 约 
束 的 优化 问题 。 由 于 在 求解 优化 问题 时 涉及 到 矩阵 求 送 ,数据 规模 大 时 计算 复杂 度 会 变 高 。 为 了 避免 大 规模 的 矩阵 求 逆 ， 
采用 拟 牛顿 法 进行 求解 ,通过 对 Hessian 矩阵 进行 分 解 ,对 计算 步 长 时 涉及 到 的 求 着 操作 进行 近似 , 且 该 算法 能 够 分 布 式 实 
项 。 仿 真实 验 表 明 ,与 现 有 算法 相 比 ,本 算法 在 大 规模 上 且 类 别 多 的 HSI 分 类 任务 下 计算 复杂 度 较 低 ,能 完成 较 高 精度 的 
1 分 类 。 
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Graph base hyperspectral images classification algorithm 
HUANG Daxin, JIANG Junzheng 
(School of Information and Communication, Guilin University of Electronic Technology, Guilin 541004, China) 


Abstract: Classification that assigns label for pixel in HSI dataset is an important pre-processed method in hyperspectral im- 
本 (CHSJI) processing，label information is useful for application such as recognition and exploration. A graph based semi- 


ervised classification method is proposed to tackle problems of large data volume, high data dimension, and small known 


ple size in HSI classification task. Dataset was modeled with graph for dimensional reducing, then the task is formulated 
.本 an unconstrained optimization problem in this method. Matrix inverse is inevitable for solving such problem, and com- 
Plexity would increase with large scale. In order to avoid large scale matrix inversion, a quasi-Newton method which approx- 
Gates inversion operation according to decomposition of Hessian matrix is used, such method can be implemented in distrib- 
uted manner. Simulations demonstrate that, compared with existing methods, proposed algorithm has lower complexity and 
higher accuracy in large scale and multi-class HSI classification task. 
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高 光谱 图 像 (hyper spectral image, 简 称 HSI) 最 
先 应 用 在 卫星 遥感 技术 上 , 现 已 在 军事 、 医 学 、 地 质 、 
农业 等 领域 得 到 广泛 的 应 用 "。 成 像 仪 对 待 测 目 标 
进行 多 窄 波 成 像 ,而 后 把 成 像 结果 进行 组 合 ,得 到 一 
个 多 频 珊 的 数据 块 , 该 数据 块 即 为 HSI。 在 HSI 的 
成 像 过 程 中 ,同一 物质 在 接收 到 不 同 波长 的 电磁 波 
时 ,会 因 内 部 材料 结构 等 因素 的 影响 而 呈现 出 不 同 的 
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反射 效果 , 相 比 于 普通 图 像 只 包含 空间 信息 ,HSI 的 
成 像 特点 使 其 携带 了 更 多 的 频谱 信息 。 对 HSI 数据 
的 研究 较为 广泛 ,有 端 元 提取 ”、 端 元 解 混 "“ ” 、 分 
类 中 等 ,其 中 分 类 是 最 为 基础 的 研究 内 容 。 分 类 的 目 
的 是 对 HSI 数据 中 每 个 像素 点 进行 类 别 标记 ,所 得 
结果 对 HSI 待 测 目标 的 识别 与 监测 等 任务 而 言 有 很 
大 的 参考 价值 "。 
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HSI 分 类 问题 可 使 用 监督 学 习 " ”算法 求解 。 
这 类 算法 在 分 类 时 需 使 用 大 量 已 标记 数据 ,然而 在 实 
际 情况 中 ,几乎 没有 原始 数据 是 已 标记 数据 ,而 对 原 
始 数据 进行 标记 需要 消耗 大 量 的 人 力 与 物力 。 因 此 ， 
能 有 效 利用 少量 已 标记 数据 与 海量 无 标记 数据 进行 
分 类 任务 的 半 监 督学 习 算法 就 受到 了 许多 关注 5 。 
而 在 分 类 问题 中 , 维度 过 高 的 特征 会 引起 * 休 斯 现 
象 "“ 维 数 灾难 ”等 问题 ,进而 影响 分 类 结果 [0 。 图 
半 监 督学 习 先 使 用 图 模型 对 数据 进行 降 维 ,接着 进行 
问题 建 模 并 优化 出 分 类 结果 ,其 中 比较 具有 代表 性 的 
有 最 小 割 算法 25 \ 局 部 全 局 一 致 算法 5 与 高 斯 域 与 
调和 函数 算法 ""。 这 些 算法 均 为 直 推 式 学 习 算 
法 5 , 相 比 于 推导 式 学 习 算法 在 优化 过 程 中 优化 的 
是 分 类 器 的 参数 53 , 直 推 式 学习 算 法 目标 函数 的 优 
化 下 接 计算 出 分 类 结果 ,因此 图 半 监 督学 习 算法 在 计 
算 速 度 上 有 一 定 有 优势 。 经 过 十 几 年 的 发 展 ,图 半 监 
入 村 习 具 有 坚实 的 理论 支撑 ,是 在 不 同 的 应 用 背景 
的 莹 用 性 都 得 到 了 验证 Ps29 。 
加 在 数据 规模 较 大 时 ,图 半 监 督学 习 算 法 的 计算 复 
架 刘 也 随 之 提升 ,针对 此 问题 ,提出 一 种 分 类 算法 ,从 
从 解 优化 问题 出 发 ,在 保证 原 有 的 分 类 精度 情况 下 降 
低 求 解 所 需 的 时 间 , 且 该 算法 可 以 分 布 式 实现 。 在 该 
外 淡 中 ,用 一 种 较为 简单 的 方式 近似 计算 目标 函数 
Hessian 矩阵 的 道 ,然后 通过 和 迭代 求 出 优化 问题 的 
名 芝 仿真 实验 表明 ,与 现 有 的 图 半 监 督学 习 算法 相 
比 $ 孩 算法 能 在 较 短 时 间 内 获得 较 好 的 分 类 结 


1 己基 础 知识 


1.1 高 光谱 图 像 简介 


与 人 眼 的 成 像 系 统 类 似 , 传 统 成 像 仪 内 部 感光 元 
器 件 只 能 感应 到 可 见 光 频段 (360 一 830 nm) 内 的 电 
磁 波 , 且 成 像 时 只 能 进行 单 宽 波段 成 像 ,因此 对 成 像 
结果 的 直观 感受 是 即 一 张 尺寸 为 m 像素 Xn 像素 
(以 下 简写 为 m Xn) 的 彩色 图 像 。 而 由 于 可 见 光 频 
段 上 的 电磁 波 均 可 以 由 红 绿 蓝 (RGB) 三 色光 组 合 表 
示 而 成 ,因此 传统 图 像 也 被 称 为 RGB 图 像 , 并 在 处 理 
时 通常 使 用 对 应 RGB 三 个 颜色 分 量 的 三 张 灰 度 图 对 
其 进行 表示 ,此 时 成 像 结 果 即 为 一 个 数据 块 ,由 3 张 
尺寸 为 m Xn 的 图 片 重 炙 而 成 站。RGB 图 像 低频 率 
分 辩 率 的 成 像 特点 ,使 得 相关 研究 极度 依赖 图 像 内 包 
含 的 空间 信息 。 因 此 对 于 RGB 图 像 而 言 ,空间 分 辨 
率 是 最 为 重要 的 特性 。 相 比 于 RGB 图 像 成 像 仪 ， 
HSI 成 像 仪 所 能 感应 到 的 电磁 波 频 段 更 为 宽广 ,覆盖 
范围 通常 为 紫外 波段 到 近 红 外 波段 。 此 外 ,HSI 成 像 


仪 另 一 重要 特性 就 是 能 进行 多 窟 波 段 成 像 , 即 成 像 仪 
在 成 像 过 程 中 ,把 待 测 物体 反射 回 的 电磁 波切 分 成 多 
个 波段 ,接着 使 用 每 个 波段 的 电磁 波 进行 成 像 , 最 后 
将 多 个 成 像 结 果 合 成 为 数据 块 。 这 样 的 成 像 方式 使 
得 HSI 以 尺寸 为 mr Xn Xd 的 数据 块 形 式 呈 现 , 其 
中 ,m Xn 是 成 像 仪 对 单个 波段 响应 的 成 像 结果 尺 
寸 ,d 则 是 成 像 时 切 分 的 波段 数 ,同一 个 像素 点 在 d 
个 频率 上 的 不 同 响应 构成 像素 点 的 频率 特性 ,频谱 特 
性 示意 图 如 图 1 所 示 。 成 像 的 特性 不 同 , 使 得 HSI 
相 比 RGB 图 像 拥有 更 为 丰富 的 频谱 信息 与 更 为 精细 
的 频率 分 辨 率 。 正 因 如 此 ,大 多 数 对 HSI 数据 的 研 
究 主要 从 其 频率 信息 出 发 。 


0.6r 
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图 1 高 光谱 图 像 频率 特性 
1.2 分 类 问题 与 建 模 


在 HSI 各 项 研究 中 ,分 类 是 较为 重要 的 研究 内 
容 。 其 根据 HSI 数据 的 频率 信息 ,对 每 个 像素 点 进 
行 类 别 标记 。 如 前 所 述 ,HSI 表示 为 一 个 尺寸 为 m xX 
n Xd 的 数据 块 ,该 数据 块 由 a 张 对 应 不 同 波段 的 图 
片 组 成 。 由 于 4 张 图 片 都 是 对 同一 个 目标 进行 成 
像 ,不 同 图 片上 的 相同 位 置 像素 点 的 值 实际 上 就 是 该 
点 对 于 不 同 频率 电磁 波 的 响应 ,而 a 个 响应 就 构成 
该 像素 点 所 对 应 物体 的 光谱 特性 。 在 HSI 分 类 问题 
中 ,光谱 特性 是 分 类 时 的 一 个 重要 参考 。 光 谱 特性 在 
引入 更 多 信息 的 同时 ,也 带 来 了 特征 的 的 维 数 问题 ， 
而 过 高 维度 会 导致 分 类 精度 下 降 。 因 此 先 使 用 图 建 
模 算法 对 其 进行 降 维 , 接 着 建立 优化 问题 ,进而 求解 
得 出 分 类 结果 。 图 是 一 种 常见 的 数据 结构 ,目前 已 在 
无 线 传感器 网 络 , 社 交 网 等 场合 中 得 到 应 用 。 图 
由 节点 与 边 组 成 ,通常 表示 为 G 二 {V,EF,W},V 与 E 
分 别 为 图 的 节点 集 与 边 集 , W 为 图 的 权 和 矩阵 , 权 和 矩 阵 
内 为 图 上 权重 值 。 节 点 之 间 通 过 边 相 连接 , 边 描述 节 
点 间 的 邻接 关系 与 相似 关系 , 边 上 权重 即 为 权 和 矩阵 中 
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的 元 素 ,两 节点 越 相 似 , 二 者 之 间 的 权重 值 越 大 。 此 
外 ,还 常用 拉 善 拉 斯 矩阵 荆 一 D 一 W 、 归 一 化 图 拉 普 
拉 斯 矩阵 二 , 一 D LD “来 表示 图 的 相关 特性 ,其 
中 DD 为 度 矩 阵 ,其 定义 为 


Ds = SW 


其 中 以 为 节点 ?i 的 邻居 集合 

一 个 像素 点 数 为 NN, 频段 数 为 4 的 HSI 数据 通 
常用 特征 矩阵 
天 一 [xiyxo wxXN] E RY 
表示 ,矩阵 中 的 每 一 行为 对 应 像素 点 的 频谱 特征 。 根 
据 该 特征 矩阵 ,用 图 对 HSI 数据 进行 建 模 ,所 建 模 图 
使 用 & 近邻 图 (kNN 图 )。 在 建 模 时 ,首先 计算 第 


(1) 


不 难看 出 ,图 半 监 督学 习 的 主要 任务 是 在 对 下, 即 
fi ,fo fe 的 求解 上 。 C 个 向 量 户 了 的 求 
解 过 程 相同 ,因此 对 仅 对 任 一 方 的 求解 进行 分 析 , 令 
了 了 一方 ， 


2 
ming(f) = 51f 了 | sf Lf 


显然 ,目标 函数 是 一 个 关于 / 数 , 且 有 人 解析 
解 , 解 为 1” = 二 (I 十 a 工 ) ”y。 但 在 数据 规模 变 大 
时 ,集中 式 求解 矩阵 的 逆 会 变 得 困难 。 


(6) 


2.1 优化 求解 
式 (6) 可 用 牛顿 法 来 求解 。 其 一 阶 与 二 阶梯 度 信 
息 分 别 为 


全 匠 点 与 余下 所 有 节点 的 欧式 距离 (这 一 距离 由 频谱 
特征 确定 ) ,接着 选择 最 相近 的 个 节点 作为 邻居 并 
与 涝 相连 ,而 后 重复 这 2 个 步 又 直到 所 有 的 节点 都 与 
其 最 近 的 k 个 邻居 相连 ,再 用 每 个 节点 与 其 的 距离 计 
BE 


A 
© 


,i EN, 
0, 其 他 情况 。 
光栅 坟 边 所 过 接 的 节 点 的 相似 程度 ,其 中 , x; 

ER™ 与 xj € RP% 分别 为 第 i 个 与 第 j 个 像素 点 的 
频 箭 特征 。 通过 以 上 步 又 , 便 得 到 一 个 kNN 图 G。 
俐 需 点 上 的 标签 用 y 了 三 Ly, ;yj 表示 , y, 代表 已 知 标 
签 6zy, 代表 未 知 标签 , 值 为 0。 在 实际 的 分 类 问题 
中 全 之 y,。 在 图 半 监 督 分 类 问题 中 ,根据 已 知 部 
分 义 签 y, 与 所 建 模 图 刻画 的 数据 特性 学 习 出 y, 的 
标签 ,完成 分 类 。 


2 分 类 问题 求解 


根据 所 建立 的 图 模型 G = {V,E,W}, 关于 HSI 
的 c 分 类 问题 可 归结 
Ding, f=E1f; — yl +L ,3) 


其 中 , y; 为 标签 one-hot 编码 矩阵 Y 的 第 j 列 , 表 示 
第 7 类 标签 在 数据 中 的 分 布 情况 .Y 的 定义 为 
1,f;=j, 
0, 其 他 情况 。 

拉 斯 矩阵 。 优 化 问题 可 以 理 
一 个 比较 符合 所 建 模 图 的 标签 one-hot 
得 到 下 后 ,通过 下 式 来 取得 最 终 的 分 


Y; = (4) 


工 , 为 图 G 的 归 一 化 拉 普 
解 为 学 习 出 一 
编码 矩阵 下 。 
类 结果 。 


(5) 


max Fi 
El 


vo(f)=a(f —y)+L,f, 
Viop(f)=al+L,. 
即 目标 函数 的 Hessian 矩阵 HH 为 xl 十 L,。 在 牛顿 
法 求解 过 程 中 ,需要 在 每 步 迭 代 中 计算 Hessian 矩阵 
的 逆 以 求 出 步 长 ;* 一 H ”YYp(f“)，, 而 这 在 矩阵 规模 
较 大 的 情况 下 会 相当 耗 时 ,在 此 提出 一 种 近似 牛 
顿 法 。 
对 于 图 归 一 化 拉 普 拉 斯 亏 ， 有 工 , 二 I 一 W,,W, 
为 图 的 归 一 化 权 和 矩阵 ,其 可 以 进一步 分 解 为 2 个 部 
分 ,对 角 部 分 Wa 与 非 对 角 部 分 Wi,， 其 中 Ws = 


(7) 


diagCW,), 且 W 一 Ww 十 W,。 基于 此 ,同样 可 将 
Hessian 和 矩阵 H 分 割 为 对 角 部 分 与 非 对 角 部 分 , 即 

H=A—B, (9) 
其 中 : 


A=al 二 (1 二 OU—W), 
B=W,, 0 —W,), 
其 中 , 0 二 0 为 控制 矩阵 分 解 的 参数 ,设置 为 1。 此 
外 ,所 建 模 的 图 都 不 带 自 环 ,因此 其 归 一 化 权 和 矩阵 的 
对 角 部 分 均 为 0, 即 Ws 一 0。 基 于 此 , 取 0=1 时 , 式 
(10) 可 进一步 更 新 为 
A 二 (2 十 oj)l， 
B—W,,+I, 
此 时 , 式 中 目标 函数 在 第 & 次 迭代 中 时 的 牛顿 步 长 为 
“一 一 (4 一 B) "yp(f*)。 根 据 文献 [21], 可 用 
s*=—(I—0QB)A vo(f’), (12) 
对 其 进行 近似 。 由 文献 L21], 0 与 Hessian 和 矩阵 的 逆 
等 效 ,并 可 以 通过 泰勒 公式 对 其 进行 近似 , 令 P = 
W, —al: 


(10) 


(11) 


0=(—P) TI 二 P， (13) 
式 (12) 中 人 为 对 角 和 矩阵 ,因此 无 需 进行 矩阵 求 道 即 
可 计算 出 步 长 ,进而 降低 了 计算 复杂 度 。 和 迭代 过 程 为 
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u’ =Be”*, 


F =A vp(f"), 
一 一 (e — Qu’), 9 
fe =f* 十 8 


2.2 分 布 式 求解 


因 所 建 图 是 ENN 图 ,基于 图 的 局 部 特性 ,上 一 
小 节 提 出 的 求解 思路 可 以 使 用 分 布 式 的 方式 来 实现 。 
在 分 布 式 计算 中 ,优化 问题 被 分 解 为 多 个 子 问题 并 分 
配 到 计算 节点 上 ,此 时 复杂 度 较 大 的 优化 问题 转换 为 
多 个 复杂 度 较 低 的 小 问题 。 分 布 式 计算 可 以 充分 有 
效 地 利用 计算 资源 , 且 很 好 地 降低 计算 复杂 度 。 分 布 
式 算法 中 ,计算 节点 除 计算 能 力 外 ,还 具有 一 定 的 通 
信 写 存储 能 力 , 且 每 个 节点 都 已 经 存储 了 图 节点 的 邻 
居 信 息 与 已 知 标签 信息 。 计 算 节点 可 以 存储 多 个 
后 就 信息 ,在 此 我 们 假设 一 个 计算 节点 仅 存储 一 个 
世 入 ,此 时 , 式 (6) 转 换 为 


也 ming(f) ek — yi): + Ff Lf ,15) 
人 > 三 


其 中 ， fj、yj; 分 别 为 1y 在 图 中 第 7 个 节点 上 的 值 。 
外 , 式 (15) 的 前 半 部 分 即 为 式 中 匹配 项 的 分 解 ， 将 

VE 加 项 分 人 为 NN 个 相互 独立 的 部 分 并 分 配 到 每 个 节 

a 分 解 后 ,节点 i 上 的 一 阶梯 度 信 息 为 

过 vp(CF) =a(fi—y)+f Zr (16) 


kb w, 为 W, 中 第 i 行 ; 列 的 元 素 , 在 图 上 则 对 应 节 
点 全 与 节点 ) 之 间 的 权重 值 。 由 此 ,节点 i 上 进行 的 
第 下 次 迁 代 则 为 

要 et =AF v9(f!), 


u’ =e? + 27et， 


j= 
9 ， (1L7) 
57 二 一 ei 十 (1 一 @)u; 4 > wt, 

jE 


油 多 


fm =fi+st。 

在 迭代 完成 后 ,所 有 的 节点 输出 的 f; 便 可 组 合 
为 one-hot 编码 矩阵 下 的 对 应 列 。 通 过 (16) 与 (17) 
可 看 出 ,节点 i 在 迭代 过 程 中 仅 需 获得 其 周边 邻居 上 
fe 以 及 的 值 即 可 计算 即 可 完成 迭代 。 对 于 节点 
i, 分 布 式 算 法 如 下 : 

while (7 过 c) do 

令 y=Y.,;， 上 k= 二 0，f"=0， 

while (| f/f — /| >E)do 

1) 与 必 中 的 节点 进行 通信 ,获取 上 了 在 必 上 的 


pa 


2) 计 算 e* 一 AP Vp(f'); 


3) 与 必 中 的 节点 进行 通信 ,获取 e 在 .上 的 值 
人 计算 


k k k 
Wr De:; 


jE 
5) 与 .WW 中 的 节点 进行 通信 ,获取 wu 在 的 值 ; 
6) 计 算 步 长 ， 


k k 
= ei 十 (1 


| k 
| Dy wt; 


jE 
DD) 计 算 Ai 二 f$ 十 5 
end while 
令 开 一 广 ; 
end while 
通过 max_F 求 得 节点 ;的 分 类 结果 了 ， 
最 终 总 合 所 有 节点 的 分 类 结果 输出 
对 算法 过 程 分 析 可 知 , 每 个 节点 在 一 次 迭代 过 程 
中 需要 与 周边 邻居 节点 进行 3 次 通信 。 


3 仿真 结果 与 分 析 


对 KennedySpaceCenter (KSC ) 与 Indian Pines 
数据 集中 进 行 分 类 仿真 。KSC 数据 内 有 202 个 波段 
图 片 ,可 分 为 13 类 ,类 别 分 布 不 平均 。Indian Pines 
内 保留 了 200 个 波段 数据 ,图 片 内 像素 分 为 16 类 ,其 
类 别 分 布 更 不 平均 ,各 类 样本 数目 差距 更 大 。 

仿真 时 ,设置 单 样 本 已 知 个 数 分 别 为 5、10、15 以 
及 20,IndianPines 的 第 7 与 第 9 类 的 已 知 数据 最 多 
为 10 个 , 且 每 种 设置 下 均 随机 生成 30 组 待 学 习 数 
据 , 最 后 将 30 组 数据 的 学 习 结 果 取 平均 作为 最 终结 
果 。 参 数 的 设置 如 表 1 所 示 。 为 了 更 好 衡量 分 类 效 
果 , 仿 真 使 用 总 体 精度 (OA)、 平 均 精度 (AA) 以 及 
Kappa 系数 三 者 的 平均 值 作为 评价 指标 。 


表 1 参数 设置 


参数 名 称 取 值 
邻居 数 5 
高 斯 核 标 准 差 c 0.1 
匹配 项 权重 = 0.5 


图 2、3 分 别 为 不 同 算法 在 不 同 数据 下 的 OA 曲 


线 图 。 从 图 2 可 看 出 ,在 KSC 数据 中 ,本 算法 的 分 类 
效果 与 文献 [19] 算 法 大 致 相同 ,高 于 文献 [20] 中 算 
法 。 从 图 3 可 看 出 ,在 IndianPines 数据 中 ,本 算法 结 
果 与 文献 [19] 中 算法 一 致 , 且 明 显 优 于 文献 L20] 算 法 
的 结 

表 2 与 表 3 分 别 为 在 各 类 已 知 样本 数量 为 15 个 
的 情况 下 ,各 算法 在 不 同 数据 下 的 运行 时 间 和 各 项 分 
类 评估 参数 。 由 于 在 不 同 已 知 样本 数 下 的 结果 分 布 
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75 | 表 3 Indian 数据 仿真 结果 
| 评价 指标 本 算法 文献 [19] 算 法 ”文献 [20] 算 法 
运行 时 间 /s 39. 10 242. 50 35.08 
03 OA/% 49.7 49. 70 8. 38 
局 60 上 AA/% 64. 41 64. 41 25. 44 
相 55 | Kappa 60. 61 60. 61 22. 36 
50 | _e 本 算法 图 4、5 分别 为 本 算法 在 不 同 数据 下 分 类 时 的 迭 
一 * 一 文献 [19] 算 法 、 外 二 ， 每 个 类 本 | 小 作 散 有 s 
天 |。 代 傅 况 。 由 于 在 每 个 类 别 进 代 的 趋势 一 致 , 任 选 其 中 
一 条 做 展示 。 通 过 2 图 可 发 现 ,本 算法 在 不 同 数据 下 
0 10 15 20 均 能 较 快 收敛 。 而 由 之 前 算法 分 析 可 知 , 本 算法 每 步 
已 知 样本 数 的 迭代 都 需要 与 邻居 进行 3 次 通信 以 完成 信息 交换 。 
J 入 、 ~ > 入 Ar、 a > £1 
2 KSC 数据 下 不 同 算法 OA 曲线 显然 ,该 算法 以 单 步 移 代 的 多 次 通信 换取 了 较 快 的 过 
代 速 度 。 
5 以 上 2 个 数据 的 仿真 结果 以 及 算法 的 分 析 可 以 
LO 60. 说 明 , 本 算法 在 数据 量 较 大 的 情况 下 更 能 体现 出 其 优 
| 势 ,更 适合 数据 量 较 大 、 类 别 数目 较 多 的 HSI 分 类 
© a 任务 。 
Dem 40T 
S 一 6 一 本 算法 16[ 
有 30 一 * 一 文献 [19] 算 法 141 
Ce 一 $ 一 文献 [20] 算 法 
CN 20| 1.2 上 
©O ,1.0F 
A Ol - 
> 1+ 所 
>< 5 10 15 20 0.6T 
到 己 知 样本 数 041| 
人 图 3 Indian Pines 数据 下 不 同 算法 OA 曲线 02| 
bs we ， ， ， 
本 加 1 2 3 4 5 6 7 
一 颈 , 仅 选取 数量 为 15 的 仿真 结果 作 展 示 。 由 表 2 渤 代 次 数 
可 看 出 ,在 KSC 数据 中 ,本 算法 比 文献 L19 算法 所 用 
。 本 图 4 ”KSC 数据 分 类 和 迭代 曲线 
时 间 大 大 减少 ,虽然 耗 时 高 于 文献 [20] 算 法 ,但 本 算 
法 各 项 评估 参数 都 表明 分 类 效果 要 更 好 。 而 表 2 中 
表明 ,在 IndianPines 数据 中 ,本 算法 的 分 类 效果 远 超 2.0 


文献 L20], 且 所 消耗 的 计算 时 间 与 文献 L20] 几 乎 一 
致 。 文 献 L19] 算 法 基于 各 类 样本 分 布 均衡 , 且 类 别 数 
目 较 少 的 前 提 , 因 此 会 在 IndianPines 数据 中 出 现 分 
类 效果 恶化 的 情况 。 


表 2 KSC 数据 仿真 结果 


评价 指标 本 算法 文献 L[19J 算 法 “文献 L[20] 算 法 
运行 时 间 /s 8. 05 35. 08 1. 35 
OA/% 73. 87 73. 87 54. 74 
AA/% 71. 13 71. 13 52. 22 
Kappa 68. 19 68. 19 48. 00 


ds 


4 
和 迭代 次 数 


6 这 


1 2 3 


5 ”Indian 数据 分 类 和 迭代 曲线 
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4 结束 语 


针对 HSI 分 类 中 未 标记 数据 多 数据 维度 高 ` 计 


算 复 杂 度 高 等 问题 ,提出 一 种 基于 图 模型 的 分 类 算 


法 , 且 该 算法 可 以 分 布 式 运行 。 该 算法 把 分 类 任务 归 
结 为 一 个 无 约束 的 优化 问题 ,首先 对 优化 问题 进行 拆 
分 ,并 将 拆 分 后 的 子 优化 问题 分 配 到 各 个 节点 上 , 接 
着 使 用 多 计算 单元 对 小 优化 问题 迭代 求解 。 仿 真实 
验 表 明 ,本 算法 在 数据 规模 大 、 类 别 多 的 分 类 任务 下 
具有 较 大 优势 。 后 续 工 作 将 使 用 图 信号 处 理 理 论 应 
用 到 数据 建 模 中 ,克服 其 对 于 样本 数目 分 布 不 均衡 的 


局 限 ; 


性 ,使 其 能 更 好 地 处 理 HSI 分 类 问题 。 
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